Realization of symmetric sharing of entanglement in semiconductor microcrystallites 

coupled by a cavity field 
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The entanglement of excitonic states in a system of spatiaUy separated semiconductor micro- 
crystallites is investigated. The interaction among the different microcrystallites is mediated by a 
single-mode cavity field. It is found that the symmetric sharing of the entanglement (measured by 
the concurrence) between any pair of the excitonic state with TV qubits defined by the number states 
(vacuum and a single-exciton states) or the coherent states (odd and even coherent states) can be 
prepared by the cavity field for this system. 
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I. INTRODUCTION 

Quantum entanglement plays an important role in the 
quantum communication and quantum-information pro- 
cessing. One can implement quantum tcleportation of a 
given state from one place to another by virtue of the en- 
tangled state 1^ . Entanglement among many particles 
is essential for most quantum communication schemes. 
The simplest generalization of the entangled states for 
more than two particles is the so-called GHZ state 
The entanglement of the nearest-neighbor pairs for an 
infinite collection of qubits arranged in a line was stud- 
ied by Wootters For a finite system, Koashi et al. ||] 
investigated the completely symmetric sharing of entan- 
glement for an arbitrary pair of TV qubits. Diir studied 
not only the symmetric sharing of the pairwise entangle- 
ment but also the nonsymmetric sharing in a system of 
TV qubits The nearest-neighbor entanglement of TV 
qubits in a ring configuration was also studied and further 
a concrete physical system of TV spin-i particles inter- 
acting via the Heisenberg antiferromagnetic Hamiltonian 
was given 0. The question arises whether it is possible 
or not to prepare states with the symmetric sharing of en- 
tanglement in some real systems within the present-day 
technology. And, if yes, how to achieve such entangled 
states? 

Any many-body system with defined qubits, if set in 
a properly chosen state, will evolve through states con- 
taining entangled qubits. Similarly, most of the ground 
states of real systems contain entangled states. However, 
for the purposes of quantum computation and quantum 
information, the most important aspects of quantum en- 
tanglement are especially (i) deterministic control over 
the quantum coherence of states and (ii) time evolution 
and occurrence of maximally entangled states. In this 



study, we focus on the latter topic, specifically, on the 
generation of the maximal pairwise entanglement. As 
was shown by Koashi et al. |^ that entanglement can- 
not be unlimitedly shared among an arbitrary number 
of qubits and the degree of bipartite entanglement de- 
creases with the increasing number of entangled pairs in 
an TV-qubit system in which any pair of particles is en- 
tangled. The maximum degree of bipartite entanglement, 
measured in the concurrence, between any pair of qubits 
is bound by 2/TV. Here, we will investigate a physical 
realization of this maximally possible bipartite entangle- 
ment. 

Within the past few years, advances in microfabrica- 
tion technology have allowed researchers to create unique 
quantum confinement, and thereby have opened up a 
new realm of fundamental physics. As low-dimensional 
semiconductor structures, quantum dots attract a con- 
siderable interest because of their atomlike properties. 
They can lead to novel optoelectronic devices that can 
be applied to the emerging fields of quantum comput- 
ing H, 1^ and quantum-information processing pT| . 
It is well known that Coulomb-correlated electron-hole 
pairs called excitons can be optically generated and con- 
trolled in a single dot , and thus can be used to store 
the quantum information and realize quantum comput- 
ing ||l^. On the other hand, a significant fraction of 
quantum computing and information schemes relies on 
the strong-coupling regime of the cavity quantum elec- 
trodynamics (QED). The observed Rabi oscillations of 
excitons in a single quantum dot |Q suggest the possi- 
bility that the quantum dot cavity QED will be realized 
in the near future. However, an essential feature of a 
quantum dot is that the electronic energy levels are com- 
pletely quantized, so the behavior of excitons deviates 
from the bosonic statistics. In the present paper, we will 
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consider some slightly bigger semiconductor microstruc- 
tures, such as the microcrystallites. In this case, the area 
of the microcrystallite is larger than that of the Bohr ra- 
dius of the exciton, and the behavior of the excitons with 
low excitation are the same as that of the bosonic parti- 
cles. Chuang et al. jl^ showed that the quantum code 
of the bosonic mode enables a more efficient error cor- 
rection. So, the mode of the excitons offers a possible 
physical implementation for such bosonic-mode coding. 
We propose a possible scheme to prepare the entangled 
excitonic states for the symmetric sharing in the system 
of N microcrystallites by virtue of the cavity QED. The 
cavity field mediates the interaction among semiconduc- 
tor microcrystallites, and then the entangled excitonic 
states can be prepared by the cavity field. 

We organize our paper as follows. In Sec. II, we will 
propose a scheme based on the present-day technology 
and model the Hamiltonian of the whole system. The 
solutions corresponding to this Hamiltonian are given 
for the general initial state. In Sees. Ill and IV, the 
bosonic exciton operator is used as an approach to deal 
with qubits uniformly. We will show how to prepare the 
entangled excitonic states with qubits defined by differ- 
ent excitonic states using various initial conditions of the 
cavity field. Any physical system cannot be isolated from 
its environment. The interaction between the system and 
the environment will result in their entanglement, then 
coherence of the qubit is destroyed with the time evolu- 
tion. So, in Sec. V, we will demonstrate the environment 
effect on the entangled states. Finally, some comments 
and conclusions will be given. 



II. MODEL AND ITS SOLUTION 

We assume that there are N spatially separated semi- 
conductor microcrystallites (also called large semicon- 
ductor quantum dots |l^) which are placed into 
an ideal semiconductor microcavity with a single-mode 
field, for example, the microcrystallites are embedded 
in a disk structure of the semiconductor, which is sim- 
ilar to the Imamoglu model for quantum dots Q. And 
we assume that the radius R of each microcrystallite is 
much larger than the Bohr radius as of excitons, but 
smaller than the wavelength A of the cavity field, that 
is, as <C i? < A. Also the distance between each pair of 
microcrystallites is larger than the optical wavelength A 
of the cavity field, and the microcrystallites indirectly in- 
teract by virtue of the cavity field. We also assume that 
there are few electrons excited from the valence band to 
the conduction band such that the exciton density for 
each microcrystallite is much smaller than the Mott den- 
sity. So, all nonlinear terms included in the interaction of 
the exciton-exciton and exciton-photon can be neglected 
in our model, and the excitons are considered as ideal 
bosons. The cavity field is assumed to resonantly inter- 
act with the zero-momentum excitons in each microcrys- 
tallite, the thermalization of the excitons is neglected. 



Under the above conditions, we can use the effective 
Hamiltonian under the rotating wave approximation as 
follows pi 
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H = hujo^a + huj^^b^-bj + fi'^gj{a'^bj + ab]), (1) 



where a(a^) is the annihilation (creation) operator of the 
cavity field with frequency oj, and bj{b^j) is the annihi- 
lation (creation) operator of the excitons in the jth mi- 
crocrystallite with the same frequency to as that of the 
cavity field. First, we assume that the coupling constants 
gj with j = 1, 2, • • • , N between the cavity field and mi- 
crocrystallites are different. We can give the Heisenberg 
equations of motion for the operators of the cavity field 
and the excitons as follows 



dAjt) 
dt 

dBj{t) 
dt 



-iY,9jBj{t), (2) 



-ig,A{t) (j = 1, 2, •• •,]¥), (3) 



where the transformations a{t) = A{t)e and bj{t) = 
Bj{t)e-''^* are applied. The solutions of Eqs. (|) and (|) 
can be obtained as 



Ait) 
Bj{t) 



a(0) cos{G't) -iJ2fjbj{0), 



(4) 



m ^ 

ifja{0), 



9j9rn[cos{G't) \ ^ , , 



G"2 



(5) 



where G" ^ ^JEU9l^ ^ /^(^) = 9jMG't)/G', 
a(0) and bj{0) {j = 1,---,N) are the initial opera- 
tors of the cavity field and excitons, respectively. We 
assume that the initial state of the whole system is 
|*(0)) = |V'(0))c|0)®^, which means that the cavity field 
is initially in the state \ip{0))c, but there is no exciton 
in any microcrystallite. Then we can obtain the wave 
function as follows 



\^{t)) = ummc\or'', 



(6) 



with the time-evolution operator U (t) 



_ -iHt/h 



III. PREPARATION OF ENTANGLED 
EXCITONIC STATE BY SINGLE-PHOTON 
STATE 

It is well known that N qubits can be defined by the 
states of N spatially separated microcrystallites. The 
two most interesting states for both experimentalists and 
theoreticians are the no-exciton and one-exciton states 
denoted by |0) and |1), respectively. So, we choose the 
computational basis states of the qubit as {|0), |1)} for 



each microcrystallite. If the cavity field is initially in 
the single-photon state \tp{0))c — a'^ |0)c, which now can 
successfully be prepared by the experiment, and no ex- 
citon is initially in any microcrystallite, then |^'(0)) = 



|0)c|0) 



Based on this initial condition, wc inter- 



polate the unit operator U^{t)U{t) into Eq. and 
consider the properties of the time-evolution operator 
UHt)OU{t) = 0{t) and U{t)\Q) = |0), the wave func- 
tion of the whole system can be obtained as follows 
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FIG. 1: Time evolutions of both the concurrence C (solid 
line) and the average photon number n — {a'' a) of the cavity 
field (dashed line) plotted for (a) iV = 3 and (b) N = 5. 



+e- 



'cos(G'i)ll)clO)^ 



(7) 



which has been returned into the original frame, and 
|l)j|0)®'^^~^'' means that — 1 microcrystallites have 
no excitons, and only one exciton is excited by the cav- 
ity field in the jth microcrystallite. We are interested 
in the entanglement between two subsystems of the ex- 
citons, such as, the nth and mth microcrystallites, then 
after tracing out the cavity field and the degrees of free- 
dom of other N — 2 microcrystallites, the reduced density 
operator for this pair of qubits can be obtained as 



Pit) 



/2jlO)(10| + /,„/„|10)(01| 
+/2|01)(01| + /„,/„|01)(10| 

{cos\G't)+ /'}|00>(00|. (8) 



-t- 



The entanglement between two qubits can mathemati- 
cally be described by using the concurrence . We as- 
sume a pair of qubits whose density matrix is pu- Then 
the concurrence of the density matrix pi2 is defined as 



C — max{Ai — A2 — A3 — A4, 0}, 



(9) 



where Ai, A2, A3, and A4, given in decreasing order, are 
the square roots of eigenvalues for the matrix 

M12 = Pl2icriy (g> (T2y)p*i2{'^ly ^ ^2y), (10) 

with the Pauli matrix 











where the asterisk denotes complex conjugation in the 
standard basis {|00), |01), |10), |11)}, and aiy and 0-23^ are 
expressed in the same basis. The entanglement of for- 
mation is a monotonically increasing function of C; and 
C = (C = 1) corresponds to an unentangled state 
(maximally entangled state). The concurrence for the 
reduced density operator (|^) can be obtained using Eqs. 
(PI) and (IlOl) as follows 



(11) 



It is found that the concurrence C periodically reaches 
its maximum value, but the values of the concurrences 
are different for different pairs, which means that the en- 
tanglements between different pairs are different. The 
coupling constants between the cavity field and micro- 
crystallites determine the entanglement of each pair. So, 
we can realize symmetric sharing of entanglement of exci- 
tonic states in semiconductor microcrystallites only when 
all microcrystallites have the same interaction with the 
cavity field, e.g., gi — §2 = ■ ■ ■ = fjN = 9, which may be 
obtained with the development of the microfabrication 
technology in the near future. Under this condition we 
can obtain the concurrence as 



C(t) = |(l-(ata)) = lsin2(Gt), 



(12) 



with G = g\l N, and the concurrence G periodically 
reaches its maximum value of 2/N. Comparing the time 
evolution of the concurrence and the average photon 
number (a^a) of the cavity field, we can easily find that 
when the average photon number is zero, the concurrence 
reaches the maximal value of 2/N for any number N of 
the microcrystallites and vice versa. As an example for 
iV = 3 and 5, Fig. |^ clearly shows this point. Under the 
condition that all microcrystallites have the same cou- 
pling with the cavity field, when the concurrence reaches 
its maximum values, the state of the microcrystallites 
system is in the generalized W state defined p^] to be 
\Wn) = iV-i/2(|10---0) -I- I010---0) + •■• +]0^--01)) 
and the cavity field is in the vacuum state, but in the 
anisotropic case, we cannot obtain the generalized W 
state for any condition. We know that a single exciton 
can be taken as a boson even for the quantum dots, so 
the assumption of the microcrystallites is not necessary 
in the case of the cavity field initially in the single-photon 
state. In the following two sections, we will mainly focus 
on the isotropic interaction of the cavity field and micro- 
crystallites to discuss the entangled coherent excitonic 
states. 
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IV. PREPARATION OF ENTANGLED 
COHERENT EXCITONIC STATE 

There are other two interesting orthogonal states called 
the even and odd coherent states (CS). These can be used 
as a robust qubit encoding for a single bosonic mode 
subject to amplitude damping, because the error caused 
by amplitude damping for this encoding can easily be 
corrected by a standard three-qubit error-correction cir- 
cuit [Q. So, in this section, we will discuss how to realize 
symmetric sharing of entanglement between any pair of 
qubits defined by the even (odd) coherent excitonic states 
in semiconductor microcrystallites. It is well known that 
one can define the even CS as the zero-qubit state |0) and 
the odd CS as the one-qubit state |1) to encode a CNOT 
quantum gate pi], that is 



|0) = N+{\a) + \-a)\ 
|1) = iV_(|a)-|-a)), 



(13) 
(14) 



with the normalization constants -/V± = (2±2e ^'"1 ) 

and I ± a) = exp(-|a|V2) X;,?=o[(±")"/^]l'^> are co- 
herent states of a bosonic annihilation operator, e.g., the 
coherent states of the annihilation operator a for the cav- 
ity field. The even and odd coherent superpositions of 
the photon states in cavity quantum electrodynamics and 
those of motional states of trapped ions can be created 
by experimentalists over the past several years. So, 
we can assume that the cavity field is initially either in 
the odd CS or in the even CS, and there are no excitons 
in any microcrystallite. In order to realize the symmetric 
sharing of entanglement, we assume that all microcrystal- 
lites have the same interaction with the cavity field, then 
the wave function of the whole system can be written by 
the factorization of the wave function [E4| as follows 



\^{t)) - iV±C/(t)[|a)c±|-a)c]|0)® 
= N±{\au{t))c\vit)a)^'' 
±\-au{t))c\-v{t)a)^''), 



(15) 



with u{t) = cos(Gi)e-*"* and v{t) = -i[sin(Gt)/ V]V]e"'"* 
and the same coupling constants between microcrystal- 
lites and the cavity field are taken. We find that all ex- 
citonic coherent states \v{t)a) in microcrystallites evolve 
periodically with time evolution and their maximal am- 
plitudes are 1/VN times the amplitude |a| of the co- 
herent cavity field. We are interested in the pairwise 
entanglement in the system of N microcrystallites. After 
tracing out the cavity field and other degrees of freedom 
for N —2 microcrystallites, the reduced density operator 
for any pair can be expressed as 

Pit) = Ni{i\vit)a){vit)a\r' + {\~vit)a){-vit)a\r^ 
±Pmv{t)a){-v{t)a\r^ 

±Pm~v{t)a){v{t)a\r'}, (16) 

where P{t) = {~u{t)a\u{t)a){{-v{t)a\v{t)a))^^^ = 
exp[-2|Q!p + 4|apsin^(Gt)/iV]. We choose the time- 
dependent even and odd CS as the basis {|0),|1)} for 



each qubit in every microcrystallite as follows |2^ 

\0)^N+mv{t)a)+\~v{t)a)), (17) 
\l)=N^mv{th)-\-v{t)a)), (18) 

where N±{t) are the normalization constants defined 
as N±{t) = (2 ± 2e-2|"l'''i"'(G0/A')-i/2. Then the re- 
duced density operator p{t) can be given, in the basis 
{|00), 1 01), |T0), |TT)}, in the following form 



Pit) = 



NUl±Pi t)), 
SNlit) 

, NUlTPjt)) 
8Nlit)Nlit) 



|00)(0 0| 4 



Nlil±Pit)) 



11) (111 



1 



NUl±Pit)) 



{|oo)(TT| + |TT)(oo|}. 



SNtit) 

-|oT)(To| + |To)(oT| 

(19) 



Following the same steps as for Eqs. (^) an d (Eo|) , we 
obtain the concurrence corresponding to Eq. (|ig|) as 



c±it) = 



g4|Q|2 sin'^{Gt)/N 

e2|"P ±1 



(20) 



where -!-(—) means that the cavity field is initially in the 
even (odd) CS. It is also found that the concurrence pe- 
riodically evolves. Although there is no simple analytical 
expression between the concurrence and the average pho- 
ton number (a'^a) of the cavity field as Eq. (p^), we find 
that if the average photon number of the cavity field 



(ata)± = |apcos2(Gi)^^|^ 



(21) 



is zero, then the value of the concurrence reaches maxi- 
mum and vice versa. So, for the fixed number N of the 
microcrystallites and the intensity |ap of the cavity field, 
the relationship of the time evolution between the con- 
currence and the average number of the cavity field is 
analogous to Fig. ^ 

It is very clear that the values of the concurrence ( [20| ) 
depend on both the number N of the microcrystallites 
and the intensity |ap of the cavity field. The maximal 
concurrence with the cavity field initially in even CS or 
odd CS decreases with the increase of the microcrystal- 
lite's number N when the intensity |ap of the cavity field 
is fixed. It is because that exp{4|ap sin2(Gi)/A^} in Eq. 
(^) is a decreasing function of the number N, so a larger 
number N corresponds to a smaller concurrence. In the 
following, we have plotted Fig. || to show the time evo- 
lution of the concurrences, and the relationship between 
the concurrence and the intensity |ap for different mi- 
crocrystallite number N. We find that the concurrence 
C-it) = C°'^'^it) periodically reaches its maximal value 
at the evolution times Gt = (2n -t- l)7r/2 (n = 0, 1, ...), 
and these points approach the upper bound value of 2/iV, 
which has been illustrated in Fig. ||(a) and Fig. ||(c) by 





FIG. 2: The concurrences are plotted as a function of time 
and intensity |ap of the cavity field when A'^ = 3 for the cav- 
ity field initially in (a) the odd CS and (b) the even CS. The 
maximal concurrences vs |ap for N=2, 3, 5, and 10 at opti- 
mum evolution times are plotted for the cavity field initially 
in (c) the odd CS and (d) the even CS. 



FIG. 3: The maximum values of the concurrences are depicted 
as a function of the number A'' of semiconductor microcrys- 
tallites for the different intensities of the cavity field with ini- 
tially the odd or even CS for (a) |Qp = 0.01, (b) |Qp = 0.1, 
(c) |Qp — 1, and (d) \af = 2, and the dashed curve corre- 
sponds to the maximum possible concurrence, given by 2/N, 
between any pair of qubits. 



an example for TV = 3, when |ap 0. It is because 
the odd CS of the cavity field is reduced to the single- 
photon state when |ap — + 0, so the concurrence C-{t) 
with initially the odd CS is reduced to Eq. ([1^ ) and 
approaches upper bound value of 2/N . Figure ||(c) also 
shows the variation for the maximal value of the concur- 
rence with the intensity |ap of the cavity field with dif- 
ferent number N of the microcrystallites, we find that the 
maximal values of the concurrence decrease with the in- 
crease of the intensity |ap, except that N = 2. However, 
when |ap 0, the even CS of the cavity field is reduced 
to the vacuum state, the concurrence C+{t) = C°™"(t) 
with initially the even CS tends to zero. The maximal 
values of C+{t) increase with the increase of the inten- 
sity lap, but when the intensity |ap is greater than a 
threshold value, which is determined by 



simplified as 



4|Q!p cosh lap 



a|2el"l' -Hcosh|a|2VF{-|a|2sech|a|2e-l"l't'^"i^l"l'} 

(22) 

for given number N {N > 2), where W{z} is the prod- 
uct log function defined as the solution for w of z = we" , 
then the concurrence gradually tends to zero. The max- 
ima of C+{t) are reached at Gt = {2n + 1)tt/2 if N and 
|ap satisfy Eq. (22). We can also find when |ap is large 



enough, such that |ap ~ 4 when > 4 and |ap w 6 
when = 3 [see Figs, ^(a) andH(b)], where bosonic ap- 
proximation for excitons is still good, then \av{t)) is ap- 
proximately orthogonal to | — av{t)) when u(t) = 0, that 
is, {—v{t)a\v{t)a) « 0. Under such condition, we can re- 
define two approximately orthogonal states |I) — \av{t)) 
and |0) — I — av{t)) as one qubit state and zero qubit 
state, then Eq. (|lj) for the reduced density operator of 
two subsystems of N {N > 2) microcrystallites can be 



p«-[|00)(00|±|ll)(ll|], 



(23) 



which means that no entanglement appears for each pair 
of microcrystallites, then the concurrences are zero for 
cavity field initially in odd CS or even CS. So when 
the cavity field is initially in the even CS, the points of 
the maximal values for the concurrence must be between 
and |ap, with condition {—v{t)a\v{t)a) « 0. These 



points are determined by Eq. (22). 

The optimum values of |ap and N maximizing C+(t) 
obtained from Eq. ( p^ ) and checked directly by numeri- 
cal maximization of Eq. (^ ) are, e.g., as follows: |ap — 
3/21n(2) « 1.04 for iV = 3, and \a\^ = ln(l + V2) ~ 0.88 
for TV — 4. While for A^ = 5 the maximum is at 
|a|2 = ln[16/9 -I- 10i/320/27 -I- 102/3l0/27]/2 w 0.81. For 
A^=6 and 8, analytical expressions can also be found. 
However, for other cases (i.e., N = 7 and N > 9) there 
are no compact-form analytical formulas for |ap corre- 
sponding to the maximal concurrence. The above numer- 
ical calculations can be seen from Figs. |(b) and |(d). 
Figure |^ also shows when the number A^ of the micro- 
crystallites is equal to 2, the maximally entangled states 
can be prepared for any intensity of the cavity field with 
initially the odd CS, but for the cavity field with initially 
the even CS, we can approximately prepare a maximally 
entangled state when the average photon number |ap is 
shghtly larger than 1, e.g., |ap > 3 [see Fig. || (d)]. 

As a comparison, we numerically show that the max- 
imal values of the concurrence in the case of the cavity 
field initially in the odd or even CS for the different num- 
ber N of the microcrystallites and the different intensity 
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lap of the cavity field from Fig. |(a) to Fig. |(d). Fig- 
ure ^ shows that the maximal values of the concurrences 
with initially the odd or even CS of the cavity field ap- 
proach each other with increasing the intensity of the 
cavity field, but they are lower than the upper bound 
values of 2/N for the concurrence. From Fig. ||(a), we 
can also find that the concurrence with the cavity field 
initially in the odd or even CS gradually approaches its 
maximal possible values of 2/N or zero when the inten- 
sity of the cavity field tends to zero, that is, |ap — > 0. 



V. DECAY OF THE ENTANGLED EXCITON 
STATES 

The quantum computation and the quantum informa- 
tion take their power from the superpositions and entan- 
glement of the quantum states, however, the necessary 
coupling of the system to the environment tends to de- 
stroy this coherence and reduces the degree of the entan- 
glement with the time evolution of the total system. So, 
in this section, we will discuss the decay of the entangled 
exciton states. For simplification of our discussion, we 
will model the interaction between the system and envi- 
ronment as follows. We assume that there is no interac- 
tion between the cavity field and environment. The dis- 
sipation of the system energy comes from the interaction 
of the excitons in the microcrystallites with environment, 
here modeled as thermal radiation fields at zero tempera- 
ture. We limit our discussion to the completely identical 
iV microcrystallites, and the coupling constants between 
microcrystallites and cavity field are the same. Under 
the above assumptions, the Hamiltonian for the system, 
environments, and their interactions can be written as 



N 



N 



N N 
j=l k j=l k 

(24) 

where a,j^k{ol f.) are annihilation (creation) operators of 
radiation fields with frequency ujj^k-, and gj_k are cou- 
pling constant between the jth microcrystallite and ra- 
diation fields. For simplicity, we assume that all gj^k are 
independent of the microcrystallite size. We assume that 
each microcrystallite separately interacts with the envi- 
ronment, but the dissipative dynamics is same for all the 
microcrystallites. The latter assumption is not necessary 
but only simplifies the degree of algebra complexity for 
calculation. 

We can obtain the Heisenberg equations of motion for 
each operator as follows 



dt 



'igA- i^gj^kAj 



fee 



(25) 



OA 
'dt 



dA 



dt 



N 



(26) 



(27) 



where the transformations a{t) = A{t)e bj{t) — 
Bj{t)e-"^\ and aj^k{t) = Aj,fc(i)e""^^-'=* are made. From 
Eq. (|^), we have 



/ dt'Bj{t')e 
Jq 



(28) 



We replace Aj_k{t) in Eqs. 
new equation as 



a) by (28), then obtain the 



—L = -igA 
dt ^ 



k 



g,.kAj,k{'^)e-'^^'---^~^' 

dt'Bj(t')e"^('^^-'=""^^*"*'H29) 



k 



We can apply the Laplace transform and the Wigner- 
Weisskopf approximation ||2^ to Eqs. (26) and (|29|), so 
that we have the solution of the cavity field as 

A{t) = u'{t)a(Q) - I J2 W&j(0) E ^jAt)^jAO), 

j J k 

(30) 

where conditions A{0) — a(0), Bj{0) = bj{0), and 
Aj,k{0) ~ CLj.kiQ) were used, and 



u'{t) 



-(7/4)* 



v'{t) = ^e-(T/4)*sin(5i)' 



6 = ./N^^^TfJW, 



sm{5t) + cos((5t)) , (31) 
(32) 
(33) 



where the small Lamb frequency shift is neglected and 
the decay rate 7 = 2TTp{LJo)\g{uJo)\'^ ■ We also use the 
former assumptions under which all microcrystallites are 
the same and have the same dissipative dynamics so that 
the decay rates 7 of each microcrystallite and the func- 
tions v'{t) of each term including operators bj{0) are the 
same. When the environment is considered, the qubits 
for each microcrystallite should be redefined as 

\0)^M+{t)[\v'{t)a) + \~v'{t)a)l (34) 
\l)=AUt)[\v'{t)a)~\~v\t)a)i (35) 

with the normalization M±(t) = [2 ± 2e-2|"'(*)"l']-i/^ 
Now, we will investigate the decay when the cavity field 
is initially in the even and odd CS, but no excitons are 
initially in any microcrystallite. After tracing out the 
degrees of the environments and other N — 2 microcrys- 
tallites for the time-dependent wave function of the whole 
system, which can also be obtained using the factorized 
form of the wave function, we get the reduced density 
operator for any two qubits as 
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FIG. 4: The time evolution of the concurrences depicted with 
the number = 3 of semiconductor microcrystalhtes for 
7/g = 0.13 (a, b) and 7/(7 = 0.5 (c, d) when the cavity field 
is initially in the odd CS (a, c) or in the even CS (b, d). 



P'±{t) 



Nl{l±P'{t)) 



8Ml{t) 
Nl{l±P'{t)) 
8Ml(t) 



00)(00| 



8Ml{t)AP{t) 



x{|01)(01| + |01)(10 
8Mlit)Ml{t) 



10)(01| + |10)(10|} 

{|oo)(TT| + |TT)(oo|}, (36) 



with P'{t) = exp [-2|a|2 (l - 2\v'{t)\^)]. Then the con- 
currences C'j.{t) corresponding to Eq. (^) can be ob- 
tained as 



p4|cnj'(t) 



1 



(37) 



where v'(t) is determined by Eq. (32). We find that the 
concurrence ( pO| ) is modified and becomes of the form 
( ^ ) after the effect of the environment is taken into ac- 
count. Now we will make further approximation. We 
assume that the couplings of the cavity field with micro- 
crystallites are stronger than the decay of the exciton, 
i.e., g 3> 7. In this case we can approximately obtain the 
concurrence Cj_(i) as 



exp 



(4|a|VA^) sin2(5i%/]V)e-(^/2)* 



1 



(38) 



Equation (^8|) shows that the entanglement between any 
pair of qubits decays in an oscillating form. It also shows 
that the increase of the microcrystallite number N results 
in the decrease of the concurrence when the decay rate 
7 and the intensity |ap are fixed. For convenience of 
discussion, we can rescale time in Eq. ( |3^ ) as t' — Gt, 
then we can find that if the coherent intensity jap of the 
cavity field and the number N of microcrystalhtes are 



given, then the larger ratio between decay rate 7 and the 
coupling constant g corresponds to the faster reduction 
of the concurrence of the entangled qubits. But if the 
ratio of 7 and g is given and the number N is fixed, 
then the higher intensity of the cavity field corresponds 
to the smaller concurrence. As an example, Fig. |j plots 
the variation of concurrence with iV = 3 for a reasonably 
good cavity j/g = 0.13 [in Fig. 4 (a) and (b)] or for 
a bad cavity 7/5 = 0.5 [in Fig. 4 (c) and (d)] according 
to Eq. (|37|). Figure || clearly demonstrates our above 
discussions. 



VI. CONCLUSIONS 

We have studied an excitonic-state implementation of 
the multiparticle entanglement based on N spatially sep- 
arated semiconductor microcrystalhtes. The interaction 
among the microcrystalhtes is mediated by a single-mode 
cavity field. We find that the entanglement (measured 
by the concurrence) between any pair of qubits that 
are defined by the excitonic number states (vacuum and 
a single-exciton states) or the coherent excitonic states 
(odd and even CS), depends on the interaction between 
the cavity field and the semiconductor microcrystalhtes. 
The entanglement between any pairs is different from one 
another for the anisotropic case. When all microcrys- 
talhtes have the same interaction with the cavity field, 
the maximal degree of the entanglement between any 
pair of qubits is the same. This condition can, prob- 
ably, be satisfied with the development of the fabrica- 
tion techniques for quantum dots and the semiconductor 
microcavity quantum electrodynamics. So, the symmet- 
ric sharing of the entanglement between any pair of N 
qubits in such a system is realizable only when the inter- 
action between N spatially separated semiconductor mi- 
crocrystalhtes and the cavity field is isotropic. Under the 
isotropic-interaction condition, when the excitonic sys- 
tem reaches maximal entanglement, all photons in the 
cavity are transformed into the excitons in the system 
of the semiconductor microcrystalhtes. The generalized 
W state and the maximal degree 2/N of entanglement 
can be obtained for the cavity field initially in the single- 
photon state. But if the cavity field is initially in the 
odd or even CS, we cannot obtain the maximal degree of 
entanglement 2 /N only except the special case where the 
cavity field is initially in the odd CS and there are two 
microcrystalhtes in the cavity |^ . 

We have also investigated the decay of any pair of the 
entangled qubits defined by the odd and even CS. The en- 
tanglement between any pair of qubits decreases because 
of the dissipation of the system energy to the environ- 
ment. If the coherent intensity of the cavity field and 
the number of microcrystalhtes are given, then with the 
rescaled time, the larger ratio between the decay rate 7 
and the coupling constant g corresponds to the faster re- 
duction of the concurrence of the entangled qubits. For 
the given ratio between the decay rate 7 and the coupling 
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constant g and the coherent intensity of the cavity field, 
the increase of the microcrystallite number N results in 
the decrease of the concurrence. But if the ratio 7/5 is 
given, and the number A'^ is fixed, then the higher inten- 
sity of the cavity field corresponds to the smaller concur- 
rence. Practically, the quality of the entanglement can 
be improved with the appearance of the new processing 
techniques and the ultrahigh finesse cavities f^. Finally, 
we should point out that our discussion is limited to the 
preparation of the entangled coherent excitonic states. 



but we cannot control them using this model. 
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